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We argue that recently proposed [Melnyk et al.. Fluid Phase Equilibr., 2009, 279, 1] a criterion to split the pair 
interaction energy into two parts, one of which is forced to be responsible the most accurate as possible for 
excluded volume energy in the system, results in expressions for the virial coefficients that improve the perfor¬ 
mance of the virial equation of state in general, and at subcritical temperatures, in particular. As an example, 
application to the Lennard-Jones-like hard-core attractive Yukawa fluid is discussed. 
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1. Introduction 

Despite the significant progress in the development of modern tools in the statistical theory of liq¬ 
uids ITlITol. there still are numerous studies where thermodynamic properties are expressed in the terms 
of virial expansions (e.g., see references 11111131 and references therein). The classical example of the 
virial expansion approach is the virial equation of state (EOS) IT^ITsl 

-^ = p + B2{T)p^ + B^(T)p^ + Bi{T)pU..., ( 1 . 1 ) 

ICbI 

where p is the pressure, /cb is the Boltzmann constant, T is the temperature, p is the number density 
and BniT), /? = 2,3,4,... are the virial coefficients. The success of the virial expansion approach relies, 
first of all, on the knowledge of the virial coefficients. The first two virial coefficients, B 2 iT) and BsiT), 
by various techniques can be obtained experimentally, while in theoretical studies both can be relatively 
easily evaluated numerically (for some fluid models even analytically (TeJ) following their definition in 
terms of Mayer cluster integrals (Ts) 

B 2 {T) = -3 j J /(ri2)dridr2 (1.2) 

and 

B3(T)^-3J j J/(ri2)/(ri3)/(r23)dridr2dr3, (1.3) 
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where 

fir)^exp[-l3uir)]-l (1.4) 

is the Mayer function and uir) is the pair interaction energy in a target model fluid. The expressions for 
higher order virial coefficients are much more complicated, especially due to a significant increase in 
the number of distinct integrals that are required to be evaluated fT7\ . Therefore, there are not so many 
potential functions uir), for which the virial coefficients BniT) of the order n>2 are known. 

Traditionally, the virial expansion approach is the most advanced for the model fluid composed of 
purely repulsive hard spheres of diameter cr. In this case, the virial coefficients Bn are independent of 
temperature and have been calculated up to the twelfth order 11811211 . Being explored in the virial EOS, 
these coefficients lead to the pressure of the hard-sphere (hs) fluid that is rather accurate in compar¬ 
ison with computer simulations data for densities up to the fluid-solid transition (22). 

However, the success of the virial expansion approach is not so evident when apart from the hard- 
sphere repulsion, the interaction potential uir) and, consequently, the exponential of the Mayer function 
in equations 111-0 both include the attractive interaction energy between molecules. In this case, the 
virial expansions approach tends to diverge when approaching the thermodynamic states associated with 
condensation. This fact imposes serious limitations on the applicability of the virial expansion approach 
to properly describe the vapour-liquid equilibrium in fluid systems. The problem of virial expansion 
divergence, in the region of condensation, turned out to be a long-standing issue in the case of Lennard- 
Jones (LJ) fluid (e.g., see recent papers by Ushcats 12311261 and references therein). 

In this communication we wish to focus on another popular model system in the liquid state theory, 
namely, the LJ-like hard-core attractive Yukawa (HCAY) fluid model (6) 


oo, 


uir) = ^ 


_^Q-zir-a) ^ 

r 


r <a, 

r^a 


(1.5) 


with za -IB [see the part (a) in figurej^for details regarding the relation between the LJ-like HCAY inter¬ 
action potential and the original LJ potential]. This model fluid has been studied intensively by computer 
experiment (2]|27]|28l as well as by other approaches, such as the mean-spherical integral equation the¬ 
ory (MSA) (6|, the MSA-based first-order perturbation theory (FMSA) (8), the MSA-based high temperature 
expansions theory I^ITol. 

As for the virial expansion approach, to the best of our knowledge, there is only one paper by Naresh 
and Singh QH, where the virial coefficients up to the sixth order, i.e., B 2 iT), B^iT ),..., B^iT), for the LJ- 
like HCAY fluid have been reported. After substituting these coefficients into the virial EOS, equation 0. 
the applicability of the latter in the case of the LJ-like HCAY fluid (see figure for illustration) can be 
summarized as follows 1131 : (i) being truncated by B^iT), the virial EOS is rather accurate in the density 
range up to a reduced density pa^ ~ 0.5, and remains to be qualitatively correct for the entire density 
range, but only for the reduced temperatures T* = /cb T/c = 1.5 and 2 that are supercritical temperatures 
for this fluid model (the critical point temperature in this case r; ~ 1.2 (7]|27]|29l); (ii) the same virial 
EOS begins to fail already right after the density pa^ ^ 0.1 in the case of the subcritical temperatures (the 
results shown by the dashed line in figure [^for the reduced temperature T* = 1). 

By using the LJ-like HCAY model fluid as a pilot system, the purpose of this exploratory study is to 
show that even in the case of the truncated virial EOS, the performance of the latter in the wide range 
of density and temperature conditions, including the subcritical ones, can be substantially improved by 
implementing the ideas that were elaborated within the framework of the augmented van der Waals 
theory (29]|30). These ideas concern the issue of a split of the total interaction potential uir) into two 
terms. Namely, in contrast to presumably van der Waals’s suggestion that the total potential energy is 
composed of the repulsion and attraction contributions, the “augmented” version of the van der Waals 
theory means that one term is representing the most accurate possible the full excluded volume energy in 
the system, that is the interaction energy between the neighbouring molecules, while the remaining part 
is responsible for the weak long-range attractive interaction energy, or the energy of cohesion, between 
the next-neighbouring molecules. 

The remainder of this paper is organized as follows: in section we provide an overview of the 
augmented van der Waals theory (29]|30). In section we discuss how the ideas of this theory can be 
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Figure 1. (Color online) Pair interaction energy in the LJ-like HCAY fluid, equation and its compari¬ 
son against the LJ counterpart (thin solid line) — part (a). The decomposition of the LJ-like HCAY pair in¬ 
teraction energy into the short-range interaction energy between two neighboring molecules or excluded 
volume interaction energy, equation with zqg = 4 [part (b)] and the weak long-range attractive in¬ 
teraction energy between target molecule and any other molecule outside the first coordination shell, 
equation {2.6) [part (c)]. 


implemented within the virial expansion approach and we present the corresponding results for the 
LJ-like HCAY fluid in section]^ We conclude with sectionj^ 


2. Augmented van der Waals theory 

Recently, it has been shown that the thermodynamics as well as the vapour-liquid equilibrium in the 
LJ-like HCAY fluid can be rather accurately described within the augmented van der Waals theory I29II30I . 
In particular, within this theory, the EOS of the LJ-like HCAY model fluid reads, 

p^po-p^a, ( 2 . 1 ) 

where, similar to the original van der Waals suggestion, coefficient a is related to the contribution from 
the attractive interaction energy between molecules, while pressure po stands for the pressure due to the 
excluded volume energy. Following van der Waals, the excluded volume pressure po originates from the 
fact that in the system of a volume V and composed of N molecules with a hard-core diameter a, each 
molecule excludes an amount of volume vq from being allowed to explore by all other molecules of the 
system. Thus, the volume accessible for molecules is reduced to V- Nvq. This phenomenon, that firstly 
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was pointed out by van der Waals, results in the pressure, 

Nk^ T pk^T 

Po =-=-> 

^ V-Nvo 1-pvo 


( 2 . 2 ) 


that is referred to as the excluded volume pressure. 

The excluded volume itself is uniquely defined upon the distance between each pair of the neigh¬ 
bouring molecules in the fiuid (see figure [^. As the first approximation, van der Waals assumed that 
the excluded volume per molecule is a constant and equals fourfold the molecular volume, i.e., vo = b = 
i2l3)7ia^. Indeed, it is the case for a dilute gaseous phase [see figure|^(b)] when the mean distance (r) 
between molecules is large (more precisely, when the mean distance (r) between the centers of the pair 
of the neighbouring molecules is larger than 2a). In the dense gaseous phase and, especially, in the liquid 
phase, the mean distance (r) between the neighbouring molecules becomes shorter than 2a, and, conse¬ 
quently, the excluded volume shells start to overlap [see figure [^(c)], resulting in the excluded volume 
per molecule being smaller than fourfold the molecular volume, i.e., vq < b. 



Figure 2. Towards definition of the excluded volume in a fluid composed of molecules with a hard core 
diameter a, and its dependence on the distance between the pair of the nearest-neighbour molecules. 


Since it is rather evident that the mean intermolecular distance (r) between the neighbouring mole¬ 
cules is affected by the number density p, it became necessary to incorporate the density dependence into 
the excluded volume vq as well. The most natural way to comply with this requirement was to utilize the 
hard-sphere fluid model for the evaluation of the excluded volume contribution to the whole spectrum of 
fluid properties and to the EOS, in particular. Such an assumption lies behind the perturbation theory of 
fluids due to Zwanzig (2) and was successfully exploited by Widom O, Barker and Henderson 0), Weeks, 
Chandler and Andersen 0 and many others |T|. 

However, although it is less evident, but the mean intermolecular distance (r) between two neigh¬ 
bouring molecules is affected by the temperature as well. Such a feature of the intermolecular distance is 
mediated by the energy of the short-range attractive interaction energy that together with the energy of 
the hard-sphere repulsion are present for the pair of neighbouring molecules. This observation suggests 
that the short-range repulsive and short-range attractive interaction energies between the neighbouring 
molecules should be incorporated into the scheme to evaluate the contribution of the excluded volume 
to the pressure. 

As the first step to comply with this idea, let us follow the suggestion (29] [30) and present the pair 
interaction energy uir) in the form, 

M(r) = Unnir) + Ma„r(r), (2.3) 

which is in contrast to the common practice (2}0 that prefers to utilize another form, 

u(r) = , (2.4) 

which assumes that the pair interaction energy uir) is separated into purely repulsive hard-sphere term 
u^^ir) and attractive interaction energy Uattrir) contributions. The interaction energy Unnir) in equa¬ 
tion (2.3) represents the full interaction energy of a target molecule and its neighbouring (nn) coun¬ 
terpart. The neighbouring molecules and the corresponding interaction energy u^nir) are identified by 
means of the range (distance) criterion. According to this criterion, the excluded volume interaction en¬ 
ergy Unnir) includes the full energy of the hard-core repulsion u^^ir) and only a part of the full attrac- 
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tion energy, namely, the part that is responsible for the interaction of a target molecule with its nearest- 
neighbour molecule only. 


Unnir) = 


U^^ir) + uf„Jr) = 


oo, 


ea 


-zoir-a) 


r < a, 


r^o. 


(2.5) 


In fact, this is the interaction energy with the molecules that belong to the first coordination shell of a 
target molecule. Following such a definition, the attraction incorporates the full attraction energy 

at the contact distance r = a between two molecules, but decays faster than the full attraction en¬ 
ergy, in order not to exceed the radii of the first coordination shell. In reality, the range of the short-range 
attraction is around one molecular hard-core diameter a, and in the present c ase can be approx¬ 
imated by fixing the decay parameter at zocr = 4. Then, the term in equation 12.31 is determined 

as the difference, uir) - Unnir), and reads 


= uir) - Unnir) = ^ 


0 , 


l^-zir-a) _^-zoir-a) 
r 1- 


r ^ C7, 

r>a. 


( 2 . 6 ) 


The pair potential u^hi^ir) corresponds to the interaction energy of the target molecule with any other 
molecule but from outside the first coordination shell. Figure shows the total pair interaction energy 
uir), the excluded volume interaction energy Unnir), 3ind the l ong-ran g e attractive interaction energy 
u^ltt^ir), all according to their definitions by equations |l.5| , l |2.5| and |2.6| , respectively. 

The shape of the long-range attractive interaction energy u^^^^ir) [see figure[^(c)] appears to be cru¬ 
cial (30) for the evaluation of the van der Waals coefficient a, that in general case is given by 


oo 

a=-2ji j g°{r)u^l^^^{r)r^dr. 
0 


(2.7) 


The function g^(r) in this equation stands now for the radial distribution function of the system with the 
excluded volume interaction potential Unnir)- We wish to stress, that only by using for excluded volume 
interaction energy Unnir) its definition accor ding to equation (2.5) , it is possible to justify the so-called 
mean-field assumption, g^(r) = 1, in equation 1 2.7 1 . In the case of the LJ-like HCAY fluid, this results in the 
simple expression. 


oo 

I = -271J I 


Jr)r^dr - -Zjia^e 


l-\-za 1 -r zod 


z^a^ 


( 2 . 8 ) 


Figure shows the results for coefficient a as they are obtained for two different choices of the ex¬ 
cluded volume model and, consequently, for two different energies of the long-range attractive interac¬ 
tion energy, the same LJ-like HCAY fluid. Namely, figure]^ (a) corre¬ 

sponds to the case when the excluded volume is described within the traditional hard-sphere model [in 
accordance with equation l |2.4) ], while figure |^(b) corresponds to the case when the excluded volume is 
described within the proposed short-range attractive Yukawa model [in accordance with equations (2.3| l 
and (2.5) 1. Two values for the coefficient a in the case of each of these models were evaluated: (i) by nu¬ 
merical integration, in accordance with the definition by equation and (ii) analytically, within the 
mean-field approximation, equation (2.8) . We note, when calculating the integral in equation (2.7) , that 
we have used for the radial distribution function g^(r) the closed-form analytical equation (3T|[32) in the 
case of hard-sphere model, and the Monte Carlo simulation data (29) in the case of short-range attractive 
Yukawa model. 

Obviously, the magnitude of the coefficient a is different for each model. However, the most intrigu¬ 
ing insight from figure comes from analysing the values of coefficient a within the same model but 
obtained from two different equations, equations (2.7) and (2.8) , respectively. We can see that in the case 
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Figure 3. (Color online) Coefficient a of the augmented van der Waals EOS (13 as calculated in ac¬ 
cordance with its definition by e quation \2.1\ (dashed line) and within the mean-field approximation, 
g^(r) = 1, given by equation (2.8) (solid line) for two different choices of the excluded volume model for 
the same LJ-like HCAY fluid: (i) excluded volume is described within the hard-sphere model — part (a), 
and (ii) within the short-range Yukawa model — part (b). 


of the hard-sphere model for excluded volume interaction [figure (a)], the mean-field and exact val¬ 
ues of coefficient a are quite different, both quantitatively and qualitatively. By contrast, in the case of 
short-range attractive Yukawa model, we could admit the tendency for coefficient a to be the same, inde¬ 
pendently of equations, \2.1) or ( |2.8) , and, consequently, which is most important — to be independent of 
the density. 


3. Augmented virial EOS 

Augmented virial EOS, which is one of the goals of the present study, will be obtained as the series in 
powers of the density p of the augmented EOS Since coefficient a does not depend on density within 
the augmented van der Waals theory, the series will concern exclusively the excluded volume pressure 
po- Namely, 


7^ = = p + BY^{T)p^ + B™(r)p3 + + ..., (3.1) 

Icb 1 Icb i 

where the virial coefficients are defined exactly as it is discussed in equations jl.2| -( |1.4| , but 

instead of the total interaction energy u{r), the short-range excluded volume interaction energy Unnir) 
must be used in the exponential of the Mayer function. 

The resulting augmented virial EOS is obtained by substituting the virial series for the excluded vol¬ 
ume pressure, equation (33> into the augmented van der Waals EOS, equation & 

^ = p + (B™(r) + ^Jp2+B™(r)p3 + B™(r)p4 + .... ( 3 . 2 ) 

From the first glance at equation \3.2) , one immediately notices two important features concerning the 
role that the long-range attraction energy plays within the augmented virial expansion approach. 

First of all, since coefficient a does not depend on the density, it follows that the long-range attraction 
energy i^am(r) contributes to the second virial coefficient only, 

(3.3) 

ICbI 

Secondly, the remaining augmented virial coefficients BsiT), B 4 iT),... all do not contain the contribution 
from the long-range attractive interaction energy being identical to those that correspond to the 
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excluded volume or short-range interaction energy Unnir) only, i.e., 

= for n>2. (3.4) 

In what follows, we apply the augmented virial EOS, equation l |3.2) , to calculate the compressibility factor 
pV/ (A/^/cb T) of the LJ-like HCAY fluid. 

4. Results and discussions 

Up to date, only the first five virial coefflcients, B^^iT), ..., B^^iT), for the excluded volume inter¬ 
action energy Unnir) are known (Ts). Figureshows a set of data for the excluded volume pressure, 
Pq = p^^, that result from the virial EOS, equation ( |3.1| , (solid lines) truncated at the sixth virial co¬ 
efficient as well as those that were obtained from computer experiment (symbols) to compare. There 
are three isotherms, namely, T* = 1,1.5 and 2 that correspond to the excluded volume pressure within 
the short-range attraction Yukawa model, while the forth isotherm represents the pressure of the hard- 
sphere fluid, i.e. corresponds to the case when T* ^ oo. The most important conclusions that follow from 
the results presented in figure [^concern the accuracy and, perhaps, even more generally — applicability 
of the virial expansions approach in the case of excluded volume interactions. First of all, we can see that 
virial EOS, equation being truncated at the sixth virial coefficient, reproduces rather accurately the 



density, 

Figure 4. (Color online) Excluded volume pressure, poKk^ T) = p^^/{k^ T), as it is modelled by the hard¬ 
core short-range Yukawa attraction energy with decay parameter zqg = A. The thick solid lines 
represent results of the virial EOS truncated at the sixth virial coefficient with coefficients B^^iT), 

..., Bq^{T) reported by Naresh and Singh ESI, while open circles are the computer experiment data by 
Shukla □ The thin dashed lines connect the symbols and are shown to guide the eye. The filled squares 
(computer experiment data) at the top and thin solid line (results of the virial EOS truncated at the tenth 
virial coefficient (18)). both correspond for the pressure, pol{k^T) = p^^l{k^T), of the hard-sphere fluid 
that represents here the high-temperature (T* oo) limit of the the excluded volume pressure, and is 
shown here for comparison purposes; the thick solid this case, like in all other cases in this figure, repre¬ 
sents the virial EOS of the hard-sphere fluid being truncated at the sixth virial coefficient. The pressure 
isotherms have been shifted for clarity. 


13501-7 






A. Trokhymchuk, R. Melnyk, I. Nezbeda 




Figure 5. (Color online) The second and third virial coefficients of the LJ-like HCAY fluid. The filled squares 
with a connecting line represent the results for BziT) and B^iT) that were obtained by Naresh and Singh 
(Is) using equations (H and (H. respectively. The open circles with a connecting line represent the 
results for B 2 {T) and B^iT) that correspond to the augmented virial expansion approach and are defined 
in accordance with equations in and {3.4} , respectively. The filled circles in part a show results for the 
second virial coefficient that are obtained in accordance with equation {3.3} but for the case when the pair 
interaction energy, u{r), is separated into two parts following a common practice EHU, i.e., it consists of 
purely repulsive hard-sphere energy, and full attractive interaction energy, Wattr(^)- 


data from computer experiment in the density range 0 < pa^ < 0.6. We note, that this observation prac¬ 
tically does not depend on the temperature; similar behaviour is found in the case of the hard-sphere 
model as well, if the virial series is truncated at the sixth virial coefficient. At the same time, by analysing 
the results of the hard-sphere model we can suggest, that truncation of the virial EOS, equation & 
at the tenth virial coefficient must be sufficient to provide rather accurate description of the excluded 
volume pressure, po = of the LJ-like HCAY fluid in the full density range. 

The definition of the augmented virial coefficients BniT) in accordance to equations {3.3} and {3.4} 
differs from a rigorous that [e.g., see equations {1.2} and {1.3} ]. Obviously, one would expect the different 
values for virial coefficients BniT) from these two definitions. Indeed, this is the case for the virial coeffi¬ 
cients B^iT), ..., B^iT), i.e., for BniT) with n>2. As an example, figure[^(b) shows the results for B^iT), 
where the augmented virial coefficient differs significantly from its conventional counterpart, remain¬ 
ing positive even for extremely low temperature. However, it does not in the case of the second virial 
coefficient when two definitions, given by equations in and {3.3} , both result in practically the same 
values of BziT) as it is illustrated in figure(a). In particular, the Boyle temperature [the temperature 
at which B 2 iT) assumes zero value] of the LJ-like HCAY fluid in both cases is the same, being fixed at ap¬ 
proximately T* ~ 2.7. We note, that this feature of the second virial coefficient is sensitive to the way how 
the total pair interaction uir) is split into the excluded volume and long-range contributions. To illustrate 
this point, figure(a) shows the results for the augmented B 2 iT) in the case when the uir) is split in 
accordance to the common practice (SJO given by equations {2.4} , i.e., when the nearest-neighbour in¬ 
teraction potential Unni^) consists of the hard-core repulsion w^r) only. We can see, that second virial 
coefficient in this case is quite different from a rigorous one. 

The results for compressibility factor, pV/iNk^ T), of the LJ-like HCAY fluid, that follow from the aug¬ 
mented virial EOS {3.2} , are shown in figure (the thick solid lines) to be compared against computer 
experiment data (7l as well as against the rigorous virial EOS O)- We can see that for all three temper¬ 
atures that include both the supercritical and subcritical conditions, there are notable improvements in 


13501-8 









8 



Figure 6. (Color online) Compressibility factor pVI{Nk^T) of the LJ-like HCAY fluid. The solid lines rep¬ 
resent the results of the truncated augmented virial EOS, equation (3.2} , truncated at the sixth virial 
coefficient, while symbols correspond to the computer experiment data by Shukla Q. The dashed lines 
show the results of the truncated conventional virial EOS, equation (H). that are taken from the study by 
Naresh and Singh (Ts). The temperature conditions are specified in the figure. We note that critical point 
temperature for the LJ-like HCAY fluid is estimated to be around T* - 1.2 (7][^[^. The curves have 
been shifted for clarity. 


searching for an agreement with computer experiment data. However, the most valuable result concerns 
the performance of the augmented virial EOS at the subcritical temperature T* = 1, where the origi¬ 
nal virial EOS equations O) fails. Some discrepancies between the augmented virial EOS and computer 
experiment data, that still are observed for densities pa^ > 0.6, are pretty similar to those that we already 
discussed in figure]^ and they should be attributed to the truncation of the excluded volume virial EOS 
(3.1) at the sixth virial coefficient. 


5. Conclusions 

In the present study, the issue of the performance of the virial expansion approach in the liquid state 
theory is discussed by using as an example the LJ-like HCAY fluid model. More precisely, we were devoted 
here to discuss the issue of the divergence of the conventional virial EOS O) at subcritical conditions 
that for the LJ-like HCAY fluid recently was reported by Naresh and Singh (13). 

To deal with this issue, the recent advances (29]|30) in the augmented version of the van der Waals 
theory have been explored. The essence of the van der Waals theory lies in the non-trivial split of total 
interaction energy into two parts which are forced to be responsible for the excluded volume energy 
and the cohesive energy, respectively. Traditionally, in what is called “van der Waals picture of liquids” 
QKD it is assumed that excluded volume part is well represented by the hard-core repulsion energy, 
while the cohesive part is associated with the remaining full attractive interaction energy. By contrast, 
within the augmented van der Waals theory (29|[30), it considers that not only the repulsive part, but the 
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full energy of interaction between the pair of neighbouring molecules, must be treated as the excluded 
volume energy, and what remains is representing the cohesion energy. 

By applying the virial expansion approach to the augmented van der Waals EOS it is obtained 
that the cohesion energy, which is the long-range part of the total interaction energy, contributes to the 
second virial coefficient, B 2 iT), only All other virial coefficients result from the short-range excluded 
volume interaction energy, Unnir), which, however, consists of both the hard-core repulsion energy and 
short-range attraction energy between the neighbouring molecules. 

To define the excluded volume interaction energy, Unnir), we have used in this study the distance 
criterion, requiring that the range of excluded volume interaction should not exceed the molecular hard¬ 
core diameter a. However, from the analysis of the second virial coefficient in figure (a) we may con¬ 
clude that as the criterion for how much of attraction energy in the total pair interaction energy, uir), 
should be included into the excluded volume term, Unnir), might be the requirement of the equality 
between the augmented second virial coefficient, equations l |3.3) , and the rigorous one, equation (TD- 
Namely, the Boyle temperature that follows from the augmented van der Waals theory should be pretty 
close to the true one, which follows from the rigorous second virial coefficient; otherwise the physics of 
two systems could differ as well. 

The range of attractive interaction energy seems to be an important issue for the convergence of the 
virial expansion approach. For example, in the limiting case of the hard-sphere repulsion, when 

there is no attractive tail at all, the virial EOS shows no sign of divergence (e.g., see reference (121 aricl 
discussion therein). Very similar conclusions can be drawn for the virial EOS l |3.11 truncated at the sixth 
virial coefficient in the case of model fiuid defined by the interaction potential Unnir). Although this 
has not been proved in the present study, we still are suggesting that virial expansions for the excluded 
volume pressure, equation (3.1) , do not diverge in the range of temperatures that are of interest for the 
parent fiuid, i.e., the LJ-like HCAY fiuid in this case, including temperatures that are below its critical point 
temperature but higher than its triple point temperature. Some discrepancies between the virial EOS and 
computer experiment data that are observed in figure|^for densities pa^ > 0.6, should be attributed to 
the truncation of the virial EOS (3.11 at the sixth virial coefficient. 

The resulting augmented virial EOS, equation (3.2) , has been tested for the LJ-like HCAY fiuid in the 
wide density range and for temperature conditions that were studied in the literature so far, including 
those where the conventional virial EOS, equation (l.l) , exhibits difficulties (13). Being rather accurate 
in the range of densities up to around pa^ = 0.6 at the supercritical temperatures, the augmented virial 
EOS, equation (3.2) , remains qualitatively correct at subcritical temperatures as well, showing no sign 
for divergence at the temperature as low as T* = 1. Nevertheless, for making final conclusion regarding 
the performance of the augmented virial EOS in the case of the LJ-like HCAY fiuid, as well as for the 
potential application of this approach to investigate more complex and/or realistic class of fiuid models, 
the evaluation of the higher order virial coefficients, namely, B^^iT ),..., Bf^iT) for the excluded 

volume interaction energy Unnir), is highly desirable. 

In general, the excluded volume pressure could be obtained by means of equation (2.2) , supposing 
that molecular excluded volume vq as function of both the density and temperature is known. Unfortu¬ 
nately, function voipy T) is not available in general case. On other hand, the pressure po can be obtained 
from the knowledge of the forces that are responsible for excluded volume. Namely, similarly to the case 
of the original van der Waals theory, when excluded volume pressure po was identified with the pressure 

of the fiuid system with a hard-sphere repulsion (r), the excluded volume pressure within the aug¬ 
mented van der Waals theory can be obtained as the pressure p^^ of the fiuid system with interaction 
potential Unnir). These data can be extracted, for instance, from computer simulation experiment (7) or 
within the integral equation theories (9][To). Such a route has been already explored (29][30|, resulting in 
the augmented van der Waals EOS for the LJ-like HCAY fiuid. The other possibility might be to utilize the 
excluded volume pressure po in the framework of the perturbed virial EOS approach I33II34I . 
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Bipia/ibHi posioiaAi^ Ta po3UJi/ipeHi/iM mctoa Ban ^ep Baayibca: 
sacTocyBaHHn pfl /1eHHapA-A>KOHciBCbKoY MOAeyii TBep^MX c((>ep 
3 npi/iT5iraHHjiM K)KaBi/i 

A. Tpoxi/iMHynpEl, p. Me;ibHi/ii^, I. He36eflcpEI 

^ iHCTHiyr (t)i3MKM KOHfleHCOBaHMX CMCTeM HAH yKpaiHM, By;i. I. CBGHLlil4bKOrO, 1, 7901 1 /IbBiB, VKpaiHa 
^ iHCTMiyr npniaiaAHoi' MareMainKM ra (t)yHflaMeHTa;ibHMx Hayx, Hai4ioHa;ibHMM yHiBepcuTei "/IbBiBCbxa 
no/iirexHiKa", 79013 /IbBie, VKpama 

^ npi/ipoAHi/iMUM c|)aKy/ibTeT, YHiBepci/iTeT iM. ^.E. flypKUHe, Vcri HaA/la6eM, 40096, MecbKa PecnyO^niKa 
^ /laOopaTopifi lepMOAHHaMiKi/i iM. E. Xa/ia, iHCTi/iTyr (t)yHAaMeHTa/ibHi/ix xiMiMHi/ix npoi^eciB AH MecbKoT 
Pecny6jiiKi/i, npara-6,16502, MecbKa PecnyO/iixa 

rioKaaaHO, u^o sanponoHOBaHi/iM He^aBHO [Melnyk et al.. Fluid Phase Equilibr., 2009, 279, 1] Kpi/irepiCi poaOi/iTT^i 
noieHi^iayiy napHoi BsaeMOAii na fl , B \ nacTi/iHi/i, o^Ha 3 jiki/ix oni/icye MoxHa TOMHiiue Bi/iK/iK)HeHi/iM oO'em b 
CMCTeMi, npMBOAWTb fl,o BMpasiB A-nfi Bipia;ibHi/ix Koe(|)ii4ieHTiB, ?iKi cyrieBo noKpaLi^yKDTb TOHHicib Bipia/ibHoro 
piBH^HH?! CTaHy B L^i/ioMy Ta TeMneparyp hi/i>kmi/ix sa Kpi/iTi/iMHy, BOKpeMa. FIk npi/iK/iaA, poBr/i^Hyro sacTO- 
cyBaHH?! fl,o /leHHapA-A>KOHciBCbKoT Mo^e/ii TBepAnx c(t)ep 3 npi/iT?iraHH?iM K)KaBi/i. 

K/iiOHOBi c/ioea: Bi/iK/ito^eHi/iM o6'€m, n/ii/iH rsepfli/ix ccpep 3 npmfiraHHfiM K)KaBi/i, Bipia/ibHe piBHBHHB crany, 
/ipyri/iM Bipia/ibHi/iM Koecl)ipi€HT, pi'bhbhhb Ban fl,ep Baa/ibca 
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